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^vq ■ Abstract. The Belle and Babar experiments have measured the branching fractions for B + — > t + v 

and D+ — » pi + V decays. From these measurements one can extract the B + and Df decay constants, 
which can be compared to lattice QCD calculations. For the Df decay constant, there is currently a 
■ 2. 1 <7 difference between the calculated value and the measured value. 
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$ : INTRODUCTION 
P-T 

Both the Belle [1] and Babar [2] experiments have measured the branching fractions 
for B + — > t + v and Df — > /i + v decays [3]. These decays proceed via the annihilation 
diagram of Fig. 1 . Within the Standard Model (SM), the predicted rates are 
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J: a{Dt^?-v) = T D ^\Vj 2 f 2 D mjm D (l-^f) . (2) 

For the Df, all parameters on the right-hand-side except for f D are well-known. The 

Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V CS \ is well-constrained by a 
global fit to several experimental observables and unitarity of the CKM matrix. Thus 
a measurement of £$(Df — > /i + v) allows one to determine the decay constant f D . This 
can be compared to QCD lattice calculations, which have become relatively precise. For 
the B + , the CKM matrix element |V> | is known to only 9%; this is nonetheless more 
precise than measurements of SS{B + — »£ + v) and allows one to extract f B . 



MEASUREMENT OF B+ -> t+v 

Belle has done two B + — > T + v analyses [4, 5]; the most recent one used 605 fb 1 of data 
and obtained evidence for a signal. This analysis employed a semileptonic tag: one B 
in an event is fully reconstructed as B -^D^°£~v, where D*°^D°K°, D°y and D° — > 
K~7t + (7t°), K~k + k~k + . The signal decays T + ^/i + v M v T , e + v e v T and T + ^7t + v T are 
then searched for by reconstructing a single track not associated with the tag side. The 
signal yield is obtained by fitting the distribution of E ECL , which is the energy sum of 
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FIGURE 1. Annihilation diagram for a heavy meson decaying to a lepton and neutrino. 



TABLE 1. Fit results for B+ -> T+ v, from Belle [5]. 
The data sample corresponds to 605 fb _1 . 
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calorimeter clusters not associated with a charged track. A peak near E ECL = indicates 
ar + -> ^ +v i^t or T+ — *■ decay. The main backgrounds are & — > c processes and 
e + e~ ^>qq continuum events. The fit is unbinned and uses a likelihood function 

* = e -^IYL»jfj( E i)> o) 

where z runs over all events (TV), j runs over all signal and background categories, 
n- is the yield of category j, and /■ is the probability density function (PDF) for 
category j. The branching fraction is calculated as SB = nj (2 • e -N B+B -), where e is 
the reconstruction efficiency as calculated from Monte Carlo (MC) simulation. The fit 
results are listed in Table 1, and the fit projections are shown in Fig. 2. The systematic 
errors are dominated by uncertainty in the background PDF and the tag reconstruction 
efficiency. The overall result is 

^+-T+v)| BeUe = (l.65l^^)xl0- 4 , (4) 

where the first error is statistical and the second is systematic. 

Babar has published two analyses of B + — > T + v decays: one using a semileptonic 
tag [6] and the other using a hadronic tag [7]. Both analyses use data samples con- 
sisting of 383 x 10 6 BB pairs. The former is similar to that used in the Belle anal- 
ysis: the tag side is reconstructed as B — > D^°£~v £ , where D*° — > D 7t ,D y and 
D° ^> K~ n + (n°) , K~k + k~iz + , K s k + k~. Babar searches for T + ^£ + v f v T , t + —>k + v x , 
and also t + — > n + 7C°v T , where for the last mode the n + n° mass is required to be near 
that of the p + . The signal is identified by plotting £e Xtra , the energy sum of calorimeter 
clusters not associated with a charged track; a peak near zero indicates T + decay. The 
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FIGURE 2. distribution of data events (points) and fit projections for B + — ► T 

(a) all t decay modes combined; (7?) t~ — >e~v e v T ; fcj t~ — > ^i~v M v T ; ft/J T~ — > 7t~ 
hatched histograms correspond to signal and background, respectively. 



V, from Belle [5]. 
V T . The open and 



signal yield is obtained by counting events in a signal region, e.g., E extm < 0.50, and 
subtracting off background as estimated from E extia sidebands. The number of events in 
the final sample is 245, the background estimate is 222 ± 13, and the resulting branching 
fraction is (0.9 ±0.6 (stat.) ±0.1 (syst.)) X 10~ 4 . 

The Babar hadronic-tag analysis is more complicated. The tagging side is recon- 
structed as B —>■ D^°X~ , where X denotes a hadronic system of total charge — 1 
composed of n^n^, n 2 (K ± ), n 3 (K®), and n 4 (n°), where n l +n 2 < 5, n 3 < 2, and 
n A < 2. The dW° is reconstructed in the same channels as those used for the semilep- 
tonic analysis, as is the T + on the signal side. A background subtraction is done on the 
tag side. The signal yield is obtained by counting events in an E extm signal region and 
subtracting off background as estimated from an E extm sideband. There are 24 signal 
candidates and 14.3 ± 3.0 estimated background events; the resulting branching fraction 
is (1.8 ±0.4 ±0.2) x 10 4 , where the first error is statistical, the second is due to 
the background uncertainty, and the third is due to other systematic sources. The data 
is shown in Fig. 3 along with projections of the fit. This result is consistent with the 



3 





extra 



0.5 

(GeV) 



FIGURE 3. E cxta distribution of data events (points) and fit projections for B + — > T + V, from Babar 
using a hadronic tag [7]. (a) t~ -^e~v e v x ; (b) T~ — > ^"V^V^ (c) T~ — ► 7C~V T ; fJJ T~ — > ^ _ 7T V T . 
The hatched histogram shows the combinatorial background component, and, for comparison, the open 
histogram shows B + — > T + v signal for a branching fraction of 0.3%. 



semileptonic-tagged result; combining the two gives 



^(fl+^T+v) = f 1.2 ±0.4 (stat.) ±0.3 (bkg.) ±0.2 (syst.)) x 1(T 4 . (5) 

Babar V / 



This is consistent with the Belle result, Eq. (4). 



MEASUREMENT OF D+ -> jU+v 

The Belle analysis of D+ — > jU + v decays [8] uses 548 fb _1 of data and searches for 
e + e~ — > DKD*n(n, y), where the primary D and ^ can be charged or neutral; the D* 
is "reconstructed" (see below) via D^y; n(n,y) =X signifies any number of additional 
pions and up to one photon (from fragmentation); the D is reconstructed via D— >Kn(n), 
where n = 1 , 2, 3 ; and neutral kaons are reconstructed via primary K is 

charged, both it and the D must have flavors opposite to that of the Df ; these constitute a 
"right-sign" (RS) sample. If the flavors are not both opposite, the event is categorized as 
"wrong-sign" (WS) and used to parameterize the background. The same classification 
applies to primary neutral K events, except for these only the D flavor must be opposite 
to that of the Df for the event to be classified as RS. 

The decay sequence is identified via a recoil mass technique. First, the recoil mass of 
the D, K, and X particles is calculated and required to be within 150 MeV/c 2 of M n *; 

then the y is included and the recoil mass is required to be within 150 MeV/c 2 of M + ; 
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FIGURE 4. Recoil mass distribution for e+e~ -^DKXyix+, from Belle [8]. 

and finally, the ji + is included and the recoil mass required to be within 0.55 GeV/c 2 
of zero. The final DKXyfi + recoil mass distribution is shown in Fig. 4; a sharp peak is 
observed near zero, indicating £>+ — > /i+v decay. 

The analysis is complicated by the fact that the recoil mass technique is very sensitive 
to the number of tracks in an event and the track reconstruction efficiency, as all tracks 
must be reconstructed for the recoil mass to be accurate. As it is difficult to simulate 
track multiplicity accurately due to uncertainties in quark fragmentation, the data is 
divided into bins of , the number of "primary particles" reconstructed in an event. 
Here, a primary particle is one that is not a daughter of any particle reconstructed in the 
event. The minimum value for is three, corresponding to e + e~ ^DKD* without any 
additional particles from fragmentation. The data is then fit in two dimensions, DKXy 
recoil mass vs. (see Fig. 5). The signal PDF is obtained from MC and modeled 
separately for different values of nj, the true number of primary particles in an event 
(rcj can differ from due to particles being lost or incorrectly assigned). 

The branching fraction is obtained from two fits: the first fit uses the DKXy recoil 
mass spectrum and yields the number of candidates; the result is N D = 32100 ± 
870 (stat) ± 1210 (syst). For this fit the background shape is taken from the WS sample 
and the background levels floated in the fit. The second fit uses the DKXy/i + recoil 
mass spectrum and yields the number of Df — > jU + v candidates; the result is iV„ v = 
169 ± 16 (stat) ± 8 (syst). For this fit the background shape is taken from a RS "D+ — > 
e + v" sample, i.e., all selection criteria are the same except that an electron candidate 
is required instead of a muon candidate. As true D+ — ► e + v decays are suppressed 
by ~ 10~ 5 , this sample provides a good model of the Df — > /i + v background. The 
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FIGURE 5. Two-dimensional DKXy recoil mass vs. n R distribution (see text), from Belle [8]. 



systematic errors listed are dominated by uncertainties in the signal and background 
PDFs and are obtained by varying the shapes of these PDFs. The branching fraction is 
the ratio N^ V /N D , corrected for the ratio of reconstruction efficiencies. The result is 

^(D+^^i+v) = (6.44 ± 0.76 (stat.) ±0.57 (syst.)) x 1(T 3 . (6) 

Belle V / 

The Babar experiment searches for D+ — > /i+v [9] using 230 fb _1 of data by 
fully reconstructing a flavor- specific dW _ , D°, or D~ decay on the tagging side. Tag 
candidates are reconstructed in the following modes: D° — > K + n~(n°) : K + K~n + n~; 
D — > K + k~k~(tiP), K®n~(n ) : K^n~n~n + , K+K n ,K°K ; D~ -> K$K~, 0p~; 
and D* ^D°7t~ with D ^^^"^ ), K$K+K , K Q s n Q . An isolated ju+ track is 
required. The neutrino momentum is taken to be the missing momentum in the event: 
p v = p e+g _ — p KSt . A photon is required and paired with the Df candidate to make a 
D* + , and the mass difference AM = M(jU + vy) -M(/i+v) is calculated. 

The data is subsequently divided into four subsamples: a tag-side mass sideband and 
a tag-side signal region for jU + and e + candidates. For both lepton samples, the tag-side- 
sideband AM spectrum is subtracted from the tag-side-signal AM spectrum (Fig. 6), and 
then the sideband- subtracted e + spectrum is subtracted from the sideband- subtracted ji + 
spectrum. The final AM distribution (Fig. 7a) is fit with signal and background PDFs; 
the signal yield obtained is N„+ v = 489 ± 55 events. 

To determine the branching fraction, the signal yield is normalized to Df — > §K + 
decays. Like the signal mode, the Df candidate is required to originate from D* + ^Dff- 
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FIGURE 6. AM spectra for the tag-side sideband region (dashed) and tag-side signal region (solid) for 
jU + (top) and e + (bottom) samples, from Babar [9]. 



The tag-side-sideband AM spectrum is subtracted from the tag-side-signal AM spectrum, 
and the resulting spectrum is fit with signal and background PDFs (Fig. 7b). The signal 
yield obtained is N^ K+ = 2093 ± 99 events. Dividing N +v by N^ n+ and correcting for 
the ratio of reconstruction efficiencies gives 



r(A + ^<^+) 



= 0.143 ±0.018 (stat.) ±0.006 (syst.). (7) 

Babar 



For this analysis, the is reconstructed via ^K + K with \M K+K _ —M^ | = AM KK < 

5.5 MeV [10]. Conveniently, CLEO has measured the branching fraction — > 
K+K-n+) for AM KK = 5 MeV; the result is (1.69 ±0.08 ±0.06)% [11]. To multiply 
the two results together to obtain SS{D^ /i + v) requires dividing Eq. (7) by — > 
K + K ) = 0.491 and subtracting (in quadrature) the 1.2% uncertainty in ^K + K~) 
from the systematic error. In addition, Babar has subtracted off a small amount of 
Df — > f Q (980)(K + K~)K + background (48 events); as this process is included in the 
CLEO measurement, these events must be added back in to Babar's §K + yield. Thus the 
Babar result becomes 



r(z>+->/i+v) 



0.285 ±0.035 (stat.) ±0.011 (syst). (8) 



&M KK =5.5 MeV 



T(D+-^K+K-7t+) 
Multiplying this by CLEO's measurement gives 

^(D+->/i+v) = (4.81 ±0.63 (stat.) ±0.25 (syst.)) x 10~ 3 . (9) 

Babar V / 
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FIGURE 7. Background-subtracted AM spectra (see text) and projections of the fit result, from 
Babar [9]. The left-most (right-most) distribution corresponds to Df —>jJ. + V (Df -^<pn + ) candidates. The 
dashed line is the background component, and the solid line is the signal plus background components 
combined. 



EXTRACTION OF DECAY CONSTANTS 

The Belle and Babar collaborations have used their measurements of 3§(B + — > x + v) and 
Eq. (1) to calculate the product of the B decay constant f B and the CKM matrix element 
\V ub \. The results are 



fa x Kb\ x 10 4 



' 9.7 ± 1.1 (stat.)!};? (syst.) GeV (Belle) 
7.2+|° (stat.) ±0.2 (syst.) GeV (Babar semileptonic) 

+2.3 



10.1 (stat.) Z l 5 (syst.) GeV (Babar hadronic). 
Taking a weighted average gives 

/ B x|^l(Belle + Babaravg.) = (9.2 ± 1 .2) X 1(T 4 GeV , (10) 
and dividing by the Particle Data Group value \V ub \ = (0.393 ±0.036)% [12] gives 

/^(BeUe+Babaravg.) = 233±37MeV. (11) 

This value is lc higher than the most recent lattice QCD results, that of the HPQCD 
collaboration (190 ± 13 MeV [13]) and that of the Fermilab/MILC collaboration (195 ± 
11 MeV [14]). 

The Heavy Flavor Averaging Group (HFAG) has calculated a world average (WA) 
value for S${Df — > /i + v) and used this to determine a WA value for the Df decay 
constant f D [15]. This value can be compared to recent lattice QCD calculations; a 
significant difference could indicate new physics. The WA value for f D is obtained by 
inverting Eq. (2): 



/* - J ^r e+V \ d2) 




G F \V cs \mAl-^)V 
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FIGURE 8. Heavy Flavor Averaging Group (HFAG) WA value for f D , from Ref. [15]. For each 
measurement, the first error listed is the total uncorrelated error, and the second error is the total correlated 
error (mostly from x D ). 



where, for &(D+ ^£ + v), the WA value is inserted. The error on f D is calculated as 
follows: values for variables on the right-hand- side of Eq. (12) are sampled from Gaus- 
sian distributions having means equal to the central values and standard deviations equal 
to their respective errors. The resulting values of f D are plotted, and the distribution is 
fit to a bifurcated Gaussian to obtain the ±1<7 errors. 

The results of this procedure are shown in Fig. 8. Also included are measurements 
of jU+v) [16] and T+v) [17] from CLEO. Thus there are three 

types of measurements: f D from the absolute D+ — > /i + v branching fraction, f D from 

the absolute £>+ — > t + v branching fraction, and f D from the — > /i + v) /T(Df — > 
K + K~n + ) ratio. The overall WA value is obtained by averaging the three results, 
carefully accounting for correlations such as the input values for \V CS \ and x D . The result 
is 256.9 ±6.8 MeV. This value is higher than the two most precise lattice QCD results, 
that of the HPQCD (241 ±3 MeV [18]) and Fermilab/MILC (249 ±11 MeV [14]) 
collaborations. The weighted average of the theory results is 241.5 ±2.9 MeV, which 
differs from the HFAG result by 2. 1 a. 
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SUMMARY 



In summary, Belle has observed B + — > T + v with 3.8c significance. From the measured 
branching fraction they determine the product f B x \ V ub \. Babar has observed B + — > T + v 
with 2.6c significance and has also measured the branching fraction to determine 
f B x |V M J. The results from the two experiments are consistent; the weighted average 
has 13% precision and is consistent with lattice QCD calculations. 

For Df — > /i + v decays, Belle has observed this mode using a recoil mass technique 
and has measured the branching fraction with 15% precision. Babar has also observed 
this mode and has measured the branching fraction relative to that for Df with 
13% precision. Dividing this by the branching fraction for (j) ^K + K and including 
Df ^> f Q (980)(K + K~)n + decays allows one to multiply by CLEO's measurement of 
3§{Df ^K + K~n + ) to obtain &(D+ — >■ ju + v). The Heavy Flavor Averaging Group has 
used the Belle and Babar measurements and also measurements from CLEO to calculate 
a world average value for f D ; the result is 256.9 ± 6.8 MeV. This value is 2.1<7 higher 
than the average of two recent lattice QCD calculations; the difference could indicate 
new physics. 
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